Theory of Initialization-Free Decoherence-Free Subspaces and Subsystems 
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We introduce a generalized theory of decoherence-free subspaces and subsystems (DFSs), which 
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large class of quantum algorithms which can function with arbitrary input states. 

PACS numbers: 03.67.Lx,03.65.Yz,03.65.Fd 



I. INTRODUCTION 

In recent years much effort has been expended to de- 
velop methods for tackling the deleterious interaction 
of controlled quantum systems with their environment. 
This effort has been motivated in large part by the need 
to overcome decoherence in quantum information pro- 
cessing tasks, affoal which was thought to be unattain- 
able at first 01 0j3j • Decoherence-free (or noiseless) sub- 
spaces 0, H ilia and subsystems 0, IE E3, E3 (DFSs) 
are among the methods which have been proposed to 
this end, and also experimentally realized in a variety of 
systems ^2 El 13 [HI ■ In this manner of passive quan- 
tum error correction, one uses symmetries in the form of 
the interaction between system and environment to find a 
"quiet corner" in the system Hilbert space not experienc- 
ing this interaction. Of the various methods of quantum 
error correction, so far only DFSs have been combined 
with quantum algorithms in the presence of decoherence 
16, 171. For a review of DFSs and a comprehensive list 
of references see Ref . . 

We have re-examined the theoretical foundation of 
DFSs and have found that the conditions for their exis- 
tence can be generalized. It is our purpose in this paper 
to present these generalized conditions. Our most sig- 
nificant result is a drastic relaxation of the initialization 
condition for DFSs: whereas it was previously believed 
that one must be able to perfectly initialize a state inside 
a DFS, here we show that this does in fact need not be so. 
Instead one can tolerate an arbitrarily large preparation 
error, which in turn means significantly relaxed experi- 
mental preparation conditions. In contrast, only a small 
preparation error can be tolerated when quantum error 
correcting codes (QECC) are used to overcome decoher- 
ence [l^ . Whether a similar generalization is possible in 
the case of QECC is an interesting open question, the 
answer to which may be within the realm of very recent 
results strengthening the DFS / QECC connection [2(| . 

The relaxation of the initialization requirement is per- 
haps most significant in light of a series of results show- 
ing that a class of important quantum algorithms (Shor 



|2ll | , Grover [23 , and Deutsch- Josza [23] included) can be 
successfully executed under im per fect initialization con- 
ditions IlllllllllllliiliilliElIll- This means 
that imperfectly initialized DFSs can be used as a "sub- 
strate" for running these algorithms. 

To present our results we first review and re-examine 
the previous results on DFSs, in Section |H] We do so 
for both general completely positive (CP) maps and for 
Markovian dynamics. The definitions we give for DFSs 
in these two cases are slightly different, reflecting the fact 
that Markovian dynamics is always continuous in time, 
whereas CP maps can also describe discrete-time evo- 
lution. In Section IIIII we present our generalized DFS 
conditions for CP maps and for Markovian dynamics. 
We illustrate the new conditions for Markovian dynam- 
ics with an example which reveals some of the new fea- 
tures. In Section IIVI we discuss the implications of our 
relaxed initialization condition in the context of quantum 
algorithms. Section [V] is devoted to a case-study of non- 
Mar kovian dynamics, intermediate between (formally ex- 
act) CP maps and (approximate) Markovian dynamics. 
A unique formulation does not exist in this case, and we 
consider the master equation introduced in Ref. 33] . The 
analytical solvability of this equation permits a rigorous 
derivation of the conditions for a DFS. For clarity of pre- 
sentation we defer most supporting calculations to the 
appendices. 



II. REVIEW OF PREVIOUS CONDITIONS FOR 
DECOHERENCE-FREE SUBSPACES AND 
SUBSYSTEMS 



We refer the reader to Ref. for a detailed review, 
including many references and historical context. Here 
we focus on aspects of direct relevance to our new results. 
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A. Decoherence-Free Subspaces 

Consider a system with Hilbert space Hs- In Refs. yj 
El IE H-l IHl a subspace "Hdfs <^ "^s was called 
decoherence-free if any state ps(0) of the system initially 
prepared in this subspace is unitarily related to the final 
state ps(t) of the system, i.e., 

ps(0) = V dPs (Q)V d Ps (t) = Up s (0)U + . (1) 

Here U is unitary and Va is the projection operator onto 
Wdfs- Important and motivating early examples of DFSs 
were given in jj, |36| [37, 38] . An alternative definition of 
a DFS is as a subspace in which the state purity is always 
one |39| : here we will not pursue this approach. 

To exploit DF-states for quantum information preser- 
vation one needs a method to experimentally verify these 
states 40], but from a theoretical standpoint one needs 
to first formulate the effect of the environment. In the 
following, we consider general CP maps and Markovian 
dynamics. 



1. Completely Positive Maps 

The modeling of environmental effects on an open 
quantum system has been a challenging problem since 
at least the 1950's ^lEHi but under certain simplify- 
ing assumptions one can obtain a simple form for the 
dynamical equations of open systems |43|. For example, 
the assumption of an initially decoupled state of system 
and bath, Psb(0) = As(0) <E> pb, results in a CP map 
known as the Kraus operator sum representation |44j |: 



Ps (t) = Tr B [A(t)(p s (0)®p B )A{tf] 
= £E a (t)p s (0)E£(t). 



(2) 



orthogonal to Hbfs by Ti DFS ±, so that lis = Hdfs © 
H UFS ± . According to Eq. (4) in j3f| the Kraus operators 
take the block-diagonal form 



E„ 



CqUdfs 






B„ 



(6) 



where the upper (lower) non-zero block acts entirely 
inside Hdfs (^dfs^); Urjps is a unitary matrix that 
is independent of the Kraus operator label a; c a is a 
scalar (J2 a \c a \ 2 = 1); and B Q is arbitrary, except that 
J2 a B^B a = I dfs _l. It is simple to verify that the DFS 
definition is satisfied in this case, with U = Udfs- 

Theorem f in [35| reads: "A subspace Hdfs is a DFS 
iff all Kraus operators have an identical unitary repre- 
sentation upon restriction to it, up to a multiplicative 
constant." This theorem is actually compatible with a 
more general form for the Kraus operators than Eq. JSJ , 
since "upon restriction to it" concerns only the upper- 
left block of E Q . We derive the most general form of E Q 
in Section IlIII below, and find that, indeed, a more gen- 
eral form than Eq. © is possible: one of the off-diagonal 
blocks need not vanish. In other words, leakage from 
7^ dfs _l into Hdfs is permitted. As we further show in 
Section ITTT1 the form © in fact appears in the context of 
unital channels. 



2. Markovian Dynamics 

The most general form of CP Markovian dynamics is 
given by the Lindblad equation 0, 0, : 



^ = -i[U s ,p s ]+C[ps], 



2 a c 



. F^ F 

2 1 ol a, 



(J) 



Here 



A(t) = T exp(-t f U(s)ds 
Jo 



(3) 



is the unitary propagator for the joint evolution of system 
and bath governed by total Hamiltonian H (T denotes 
time-ordering and we work in units such that K = 1); the 
"Kraus operators" {E Q } are given by 



E a = VA„(//|A 



(4) 



where are bath states in the spectral decompo- 

sition pb = J2 V Aj,|f)(z/|. Trace preservation of ps(t) 
implies the sum rule 



E E " E ° 



(5) 



where F„ are bounded (or unbounded, if subject to ap- 
propriate domain restrictions 0> HH) operators ac ting 
on 7^s, and where Hg may include a Lamb shift [50|. 
Given such dynamics, one restores unitarity [i.e., the DFS 
definition with U generated by the Hamiltonian H5] 
if the Lindblad term £[ps] can be eliminated. According 
to Refs. [H EDI j a necessary and sufficient condition for 
this to be the case is 

F a |t) = c a \i), (8) 

where Hbfs = Span{|i)} and {c a } are arbitrary complex 
scalars. Thus the Lindblad operators can be written in 
block-form as follows: 



Cq<I A a 
B a 



(9) 



where la is the identity operator on the system. 

In |35| a DFS-condition was derived for general CP 
maps of this type. We denote the subspace of states 



with the blocks on the diagonal corresponding once again 
to operators restricted to Hufs and 7i DFS i . Note the ap- 
pearance of the off-diagonal block A Q mixing 7Ydfs and 
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H dfs j_; its presence is permitted since the DFS condi- 
tion (JSJ gives no information about matrix elements of 
the form (i\F a \j ± ), with \i) g Hdfs and Ij- 1 ) G 7i DFS i. 

As observed in Refs. JEHU, one should in addition re- 
quire that Hs does not mix DF states with non-DF ones, 
i.e., mixed matrix elements of the type (j'^IHsli), with 
|i) € Hbfs and Ij- 1 ) € W DFS -l, should vanish. We show 
below that this condition must be made more stringent. 



B. Noiseless Subsystems 

An important observation made in Ref. [E] is that there 
is no need to restrict the decoherence-free dynamics to a 
subspace. A more general situation is when the DF dy- 
namics is a "subsystem" , or a factor in a tensor product 
decomposition of subspace. Following Ref. Jg| , this comes 
about as follows. Consider the dynamics of a system S 
coupled to a bath B via the Hamiltonian 

H = H s ®I B +I s <g)H B + H / , (10) 

where H5 (H B ), the system (bath) Hamiltonian, acts 
on the system (bath) Hilbert space Hs (Wb); Is (Is) is 
the identity operator on the system (bath) Hilbert space; 
Hi is the interaction term of Hamiltonian which can be 
written in general as J2 a ^ a ® ^ the system Hamil- 
tonian Hs and the system components of the interaction 
Hamiltonian, the S a 's, form an algebra S, it must be f- 
closed to preserve the unitarity of system-bath dynamics. 
Now, if A is a f-closed operator algebra which includes 
the identity operator, then a fundamental theorem of C* 
algebras states that A is a reducible subalgebra of the full 
algebra of operators [S^ ■ This theorem implies that the 
algebra is isomorphic to a direct sum of dj x dj complex 
matrix algebras, each with multiplicity n,j: 

5=0W®M(dj,C) (11) 
JeJ 

Here J is a finite set labeling the irreducible components 
of <S, and A4(dj, C) denotes a dj x dj complex matrix al- 
gebra. Associated with this decomposition of the algebra 
S is a decomposition of the system Hilbert space: 

Hs=0C" J ®C' lj . (12) 
JeJ 

If we encode quantum information into a subsystem 
(factor) C" J it is preserved, since the noise algebra S 
acts trivially (as I nj ). In such a case C nj is called a 
decoherence-free, or noiseless subsystem (NS) Q. Ex- 
amples of this construction were given independently in 
Refs. HEt 



1. Completely Positive Maps 

As the Kraus operators are given by Eq. Q , they take 
the form of the decomposition (lll|l : 

E a = I nj ® M a (dj), (13) 
JeJ 

where M a (dj) is an arbitrary <ij-dimensional complex 
matrix. Therefore a factor C™ J is a NS if the Kraus 
operators have the representation ljT3|) . 



2. Markovian Dynamics 

The aforementioned reducibility theorem [^2 does not 
apply directly in the Markovian case, since the set of 
Lindblad operators {F Q } need not be closed under con- 
jugation. Nevertheless, as shown in |l(j, the concept of 
a subsystem applies in the Markovian case as well: the 
condition for a NS was found to be 

F a V d = Inj <g> M a (dj)V d , (14) 

with the M Q again being arbitrary complex matrices 
and Vd being the projection operator onto a given sub- 
space C nj (g) C dj . The NS is then a factor C" J as in 
Eq. (O, with the same tensor product structure as in 
Eq. O. 



III. GENERALIZED CONDITIONS FOR 
DECOHERENCE-FREE SUBSPACES AND 
SUBSYSTEMS 

We now proceed to re-examine the conditions for 
the existence of decoherence-free subspaces and sub- 
systems. We will show that the conditions presented 
in the papers l aying th e general theoretical foundation 
[E IE C3 I35L l5l| , can be generalized and sharp- 
ened, both for CP maps and for Markovian dynamics. 
Our main new finding is that the preparation step can 
tolerate arbitrarily large errors. Relatedly, we consider 
the possibility of leakage from outside of the protected 
subspace/subsystem into it. Previous studies did not al- 
low for this possibility, but we will show that it can be 
permitted under appropriate restrictions. In doing so we 
generalize the definition of a NS with respect to the orig- 
inal definition that relied on the algebraic isomorphism 
(|llf) (see Ref. 20] for a related recent result). In the case 
of Markovian dynamics, our main new finding is that if 
one demands perfect initialization into a DFS then the 
condition on the Hamiltonian component of the evolution 
is modified compared to previous studies. 

The derivation of these results is somewhat tedious. 
Hence, for clarity of presentation we focus on presenting 
our generalized conditions in this section. Mathematical 
proofs are deferred to the appendices. We begin with 
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the simpler case of decoherence-free subspaces and con- 
sider the case of CP maps and Markovian dynamics. We 
then move on to the case of decoherence-free (noiseless) 
subsystems. The case of non-Markovian continuous-time 
dynamics is treated later, in Section Ivl 

A. Decoherence-Free Subspaces 

The system density matrix ps is an operator on the 
entire system Hilbert space Tls, which we assume to be 
decomposable into a direct sum as H = Ti.F>FS ( S'hir)FS ± ■ 
It is convenient for our purposes to represent the system 
state (and later on the Kraus and Lindblad operators) in 
a matrix form whose block structure corresponds to this 
decomposition of the Hilbert space. Thus the system 
density matrix takes the form 

/ PDFS 92 \ n 

We also define a projector 

Vdfs = ( Idfs ) , (16) 
so that p DFS = T'bfsPs'Pbfs- Finally, 

are projection operators onto Hbfs and Tl DFS ±- , respec- 
tively. 

1. Completely Positive Maps 

The original concept of a DFS, Eq. poses a practi- 
cal problem: the perfect initialization of a quantum sys- 
tem inside a DFS might be challenging in many cases. 
Therefore we introduce a generalized definition to relax 
this constraint: 

Definition 1 Let the system Hilbert space Tis decom- 
pose into a direct sum as TL = TiDFsffiT^DFS^ > an d P ar ~ 
tition the system state ps accordingly into blocks, as in 
Eq. Assume Pdfs(O) = 7 ? dfsAs'(0)'Pd FS / 0. Then 

Wdfs * s called decoherence-free iff the initial and final 
DFS-blocks of ps are unitarily related: 

Pdfs (t)=U D FSPDFs(0)Uj )FS) (18) 

where Udfs * s a unitary matrix acting on Hdfs • 

Definition 2 Perfect initialization (DF subspaces): 
P2 = and P3 — in Eq. y3J). 

Definition 3 Imperfect initialization (DF subspaces): 
pi and/ or p^ in Eq. k!5\) are non-vanishing. 

We prove in Appendix IA II 



Theorem 1 Assume imperfect initialization. Let U be 
unitary, c a scalars satisfying \c a \ 2 = 1, and B a arbi- 
trary operators on 7i DFS i satisfying ^ Q BjjB a = Irj F gj- . 
A necessary and sufficient condition for the existence of a 
DFS with respect to CP maps is that the Kraus operators 
have a matrix representation of the form 

e - = ( c o U b°J- < M > 

This form is identical to the previous result @, with 
the important distinction that due to the new definition 
of a DFS, Eq. i|18|) . the theorem holds not just for states 
initialized perfectly into 7Ydfs, but for arbitrary initial 
states. Note that unlike fault-tolerant QECC, where the 
initial state must be sufficiently close to a valid code state 
|19| . here the initial state can be arbitrarily far from a 
DFS-code state, as long as the initial projection into the 
DFS is non-vanishing. 

These observations lead us to reconsider the original 
definition, wherein the system is initialized inside the 
DFS. This situation admits more general Kraus opera- 
tors. Specifically, we prove Appendix IA ll that: 

Corollary 1 Assume perfect initialization. Then the 
DFS condition is: 

E a = ( C f £ ) , (20) 

where U is unitary. 

Note that due to the sum rule ^2 a EjjE Q = I the 
otherwise arbitrary operators A a and B Q satisfy the 
constraints (i) J2 a ^ a A a + B^Bq = I DFS ± and (ii) 
^2 a c* A„ = 0, and where additionally the scalars c a 
satisfy (iii) Y, a l c «| 2 = !• 

In contrast to the diagonal form in the previous condi- 
tions © and l|19fl , Eq. (|20() allows for the existence of the 
off-diagonal term A a , which permits leakage from Ti. DFS ± 
into Hufs- This more general form of the Kraus oper- 
ators imply that a larger class of noise processes allow 
for the existence of DFSs, as compared to the previous 
condition (JSJ). 1 



2. Unital Maps 

A unital (sometimes called bi-stochastic) channel is a 
CP map 3>(p) = ^] a E Q pE] 1 that preserves the identity 
operator: €>(I) = J^EaE^ = I. Consider the fixed 
points of 3>, i.e., Fix(<I>) = {p : 3?(p) = p}. Such states, 



We re-emphasize that Theorem 1 in l.'-iHl is compatible with 
Eq. I2UI : the latter generalizes the explicit matrix representa- 
tion Eq. (4) given in that paper [condition EJ in the present 
paper], but does not invalidate Theorem 1 in l35l . 
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which are invariant under are clearly examples of DF- 
states of the corresponding channel. 

Recently it has been shown that the fixed point set of 
unital CP maps is the commutant of the algebra gener- 
ated by Kraus operators [53j . In other words, if £ is the 
set of all polynomials in {E Q }, or £ = Alg{E Q }, then 



Fix(#) = {TeS(H):[T,£]=0}, 



(21) 



where B(H) is the (Banach) space of all bounded oper- 
ators on the Hilbert space Tt. In other words, the fixed 
points of a unital CP map, which are DF states, can alter- 
natively be characterized as the commutant of Alg{E Q }, 
i.e., the set {T}. It is our purpose in this subsection to 
show that, under our generalized definition of DFSs, this 
characterization of DF states is sufficient but not neces- 
sary. 

Consider the generalized DFS-condition lt2*U|) applied 
to unital maps. We have 



*M = E( cJ o dfs bIM C ^ FS b°.)- < 22 > 

Unitality, <&(I) = I, together with J2 a \°a\ 2 = 1 implies: 

( Idfs + E„ A Q A[ Ea \ _ T m] 

This implies the vanishing of the matrices A Q , so that 
we are left with the Kraus operators in the simple block- 
diagonal form: 



E„ 



c a I 
B Q 



(24) 



together with the additional constraint J2 a ^aB„ = 
I dfs j_ (which, in the present unital case, naturally sup- 
plements the previously derived normalization constraint 
^2 a B^B Q = I dfs j_). Thus, unitality restricts the class 
of Kraus operators, so that in fact we must assume the 
DFS-condition rather than (5Djl. This then means 
that we may consider the generalized DFS definition 
Eq. P|) . 

Next, let us find the commutant of this class of Kraus 
operators. First, 



Alg{E Q } = { 



poly(c Q )I 







poly(B Q ) 



}• 



(25) 



where poly(a;) denotes all possible polynomials in x. Rep- 
resenting an arbitrary operator T S B(H) in the form 



L M 
N P 



(26) 



it is simple to derive that the commutant of Alg{E Q } is 
the space of matrices T of the form 



L 

cl 



(27) 



where L and c are arbitrary. The aforementioned the- 
orem 153] states that the fixed-point set of the channel, 
i.e., the DF states, coincides with this commutant. Of 
course, for T to be a proper quantum state it must be 
Hermitian and have unit trace, whence c > and L is 
Hermitian. Subject to these constraints we see that the 
aforementioned theorem [53j gives a sufficient, but not 
necessary characterization of the allowed DF states. In- 
deed, the form l|2Tjl arises as a special case of our consid- 
erations, where we allow for T to be a state with support 
in 7i DFS j_ , but not of the most general form allowed by 
Eq. jTSjl , which includes off-diagonal blocks. 



3. Markovian Dynamics 

In the case of CP maps we are only interested in the 
output state and the intermediate-time states are ig- 
nored. Since, as is well known, Markovian dynamics is 
a special case of CP maps (e.g., 0, EH)> one mav 01 
course apply the results we have obtained above for gen- 
eral CP maps in the Markovian case as well, provided one 
is only interested in the state at the end of the Markovian 
channel. However, one may instead be interested in a 
different notion of decoherence-freeness, wherein the sys- 
tem remains DF throughout the entire evolution. Such 
a notion is more suited to experiments in which the final 
time is not a priori known. This is the notion we will pur- 
sue here in our treatment of continuous-time dynamics, 
in both the Markovian and non-Markovian cases. Thus, 
while we allow that the system not be fully initialized 
into the DFS, we require that the component that is, un- 
dergoes unitary dynamics at all times. Correspondingly, 
we define a DFS in the Markovian case as follows: 

Definition 4 Let the system Hilbert space Tls decom- 
pose into a direct sum as Tls — ^DFsffi^DFS^ > an d P ar ~ 
tition the system state ps accordingly into blocks. Let 
Vufs be a projector onto 7Ydfs and assume /9dfs(0) = 
'PdfsPs(0)'Pbfs 0- Then Hufs *s called decoherence- 
free iff pdfs undergoes Schrddinger-like dynamics, 



dp 



DFS 



dt 



= — «[H D fs, PdfsL 



(28) 



where Hdfs * s a Hermitian operator. 



Before presenting the DFS conditions, let us recall 
the quantum trajectories interpretation of Markovian dy- 
namics [3 EM EH ■ Expanding Eq. Q) to first order in 
the short time-interval r yields the CP map 



(3=0 



where 



W - I- i rH s --^Ft e F Q , 

a 

W^ >0 = VrFf}, 



(29) 

(30) 
(31) 
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and to the same order we also have the normalization 
condition 



E 

0=0 



wlw fj = I. 



(32) 



Thus the Lindblad equation has been recast as a Kraus 
operator sum , but only to first order in r, the coarse- 
graining time scale for which the Markovian approxi- 
mation is valid [Bol ]. This implies a measurement in- 
terpretation, wherein the system state is ps(t + t) = 
Wpp(t)Wt/pp (to first-order in r) with probability pp — 

Tr[Wpp(t)Wp]. This happens because the bath func- 
tions as a probe coupled to the system while being sub- 
jected to a quasi-continuous series of measurements at 
each infinitesimal time interval r l33|. The result is the 
well-known quantum jump process0 EH EH > wherein 
the measurement operators are Wo « exp(— irH c ), the 
"conditional" evolution, generated by the non-Hermitian 
"Hamiltonian" 



Hs-l^FtFc 



(33) 



and t/rFp (the "jump"). Note that Hg is here meant 
to include all renormalization effects due to the system- 
bath interaction, e.g., a possible Lamb shift (see, e.g., 
Ref. E3|)- By a simple algebraic rearrangement one can 
rewrite the Lindblad equation in the following form: 



Ps 



= -i(H cPs - psrHt) + F aPs Fl 



(34) 



where according to the above interpretation the first term 
generates non-unitary dynamics, while the second is re- 
sponsible for the quantum jumps. 

Now recall the Markovian DFS condition derived in 
Refs. [H E^: the Lindblad operators should have triv- 
ial action on DF-states, as in Eq. JSJ), i.e., F a \i) — 
c a \i). Viewed from the perspective of the quantum-jump 
picture of Markovian dynamics, this implies that the 
jump operators do not alter a DF-state, i.e., the term 
Y, a F <*PS F l m Eq- transforms p s to £ Q \c a \ 2 p s 

and thus has trivial action. 

Given Eq. JSJ, the Lindblad operators can be written 
in block-form as follows [Eq. J3J]: 



c a I -A_ a 
B„ 



(35) 



with the blocks on the diagonal corresponding once again 
to operators restricted to Hdfs and 7i DFS i . Note the ap- 
pearance of the off-diagonal block A Q mixing Hdfs and 
7i DFg j_ ; its presence is permitted since the DFS condi- 
tion (jSJ gives no information about matrix elements of 
the form (zIFqIj- 1 ), with \i) G Hdfs and \j ± ) G H DFS ±. 

As observed in [6| , one should in addition require that 
Hs does not mix DF states with non-DF ones. It turns 
out that this condition is compatible with the case that 



the DF state is imperfectly initialized (Definition 0) . In 
this shown in Appendix IA 21 the following theo- 

rem holds: 

Theorem 2 Assume imperfect initialization. Then 
a subspace Hdfs of the total Hilbert space TL is 
decoherence-free with respect to Markovian dynamics iff 
the Lindblad operators F„ and the system Hamiltonian 
Hs assume the block- diagonal form 



H. 



Hdfs 






Hdfs j 



c a I 
B„ 



where Hdfs o,nd H DFS i are Hermitian, c 
and B Q are arbitrary operators on 7i DFS J_ . 



(36) 



scalar s, 



But, as is clear from the quantum jumps picture, 
in particular Eqs. (|33|l . (|34(l . there also exists a non- 
Hermitian term, which appears not to be addressed prop- 
erly by merely restricting H5. Indeed, this is the case if 
one demands that the system state is perfectly initialized 
into the DFS (Definition |5J). As shown in Appendix I A 21 
the full condition on the Hamiltonian term then is: 



■iH s + i^FtF a )|fc x ) =0, 



Vz,/c J 



(37) 



where |i) G Hdfs, \k x ) G H UFS ±. Applying the DFS 
conditions (|9"|). H37|) . the Lindblad equation J7J) reduces 
to the Schrodinger-like equation Ij28(l . Combining these 
results, we have: 

Theorem 3 Assume perfect initialization. Then a sub- 
space Hdfs °f the total Hilbert space TL is decoherence- 
free with respect to Markovian dynamics iff the Lindblad 
operators F Q and Hamiltonian H5 satisfy 



B„ 



'PdfsHs'P 



DFS 



c a A n . 



(38) 
(39) 



Note that Hg (which, again, includes the Lamb shift) 
must satisfy a more stringent constraint than previously 
noted due to the extra condition on its off-diagonal block. 
This has implications in examples of practical interest, as 
we next illustrate. 



4- Example (significance of the new condition on the 
off-diagonal blocks of Hs ) 

We present an example meant to demonstrate how the 
new constraint, Eq. (|37|) [or, equivalently, Eq. H39JI ]. may 
lead to a different prediction than the old constraint, that 
matrix elements of the type (j |Hs|i), with \i) G Hdfs 



and G H 



DFS A 



should vanish. 



Consider a system of three qubits interacting with a 
common bath. The system is under influence of the bath 
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via: 1) Spontaneous emission from the highest level |111) 
to the lower levels, 2) Dephasing of the first and the sec- 
ond qubits. For simplicity we set the system and bath 
Hamiltonians, Us and Hb, to zero. The total Hamilto- 
nian then contains only the system-bath interaction: 

H/ = +<7f)®B + A 2 [(0-7/+o7r+o-.r)(g)b t 

+ +<t+ + a+)®b], (40) 

where 

erf = |001) (111 |, erf = |010)<111], (Tf = |100>(111|. 

(41) 

and b is a bosonic annihilation operator. 

The corresponding Lindblad equation may be derived, 
e.g., using the method developed in Ref. [5(j. It may 
then be shown that 

1 2 

C \Ps\ = 2 E[ F "^ F !] + FiPsMl (42) 
i=i 

where the Lindblad operators are 

Fi = \/dk(uuKi + U12K2), 

F 2 = v^2(w2iKi +m 22 K 2 ). (43) 

Here Ki = af + <j|, K 2 = erf + erf + <rf , and {d\, d 2 } 
are the eigenvalues of the Hermitian matrix A = [a,j] 
of coefficients in the pre-diagonalizcd Lindblad equation, 
with the diagonalizing matrix denoted U = [ity]. 

Now let us find the DFS conditions under the assump- 
tion of perfect initialization. The previously-derived 
Eq. © yields that {|000), |001)} is a DFS, since K 2 an- 
nihilates these states, and they are both eigenstates of 
Ki with an eigenvalue of +2: 

Fi|000) = 2VdiMii|000), F 2 |000) = 2-^^21 1 000) 
Fi|001) = 2v/rfiuii|001), F 2 |001) = 2VdiUn|001). 

(44) 

However, the new condition l|37|) tightens the situation. 
Choosing as representatives the states 1 001) G 7Ydfs and 
I 111) S H DFS ±, we find from Eq. JSJ: 

2 

(001|^FtF a |lll) - 2d 1 u* n u 12 + 2d 2 u* lU 22 

= 0. (45) 

Since u^un + it2 1 u 22 = (from unitarity of U), we see 
that the new condition imposes the extra symmetry con- 
straint di = e? 2 . This example illustrate the importance 
of the new condition, Eq. I|37|) . 

B. Noiseless Subsystems 

We now consider again the more general setting of sub- 
systems, rather than subspaces. 



1. Completely Positive Maps 

Suppose the system Hilbert space can be decomposed 
as TLs = Hns <8 Tti n © W ou t , where Hns is the factor in 
which quantum information will be stored. The subspace 
H ou t may itself have a tensor product structure, i.e., ad- 
ditional factors similar to 7Yns may be contained in it 
[as in Eq. (|12(l ]. but we shall not be interested in those 
other factors since the direct sum structure implies that 
different noiseless factors cannot be used simultaneously 
in a coherent manner. As in the DF subspace case con- 
sidered above, we allow for the most general situation of 
a system that is not necessarily initially DF. To make 
this notion precise, let us generalize the definitions of the 
projector "Pdfs and projection operators VdjV^ given 
in the DFS follows: 

^NS-in = ( Ins ® I in ) , (46) 
/ Ins <8> Iin \ _ ( \ 

Vd -{ 0)' v ^-{o i NS ®i in J ( 4? ) 

There is no risk of confusion in using the DFS notation, 
Pd, for the NS the DFS case is obtained when 

Iin is a scalar. 

The system density matrix takes the corresponding 
block form 

Definition 5 Let the system Hilbert space 7is decompose 
as TCs = Hns ® Ti-in ffi Hout, and partition the system 
state ps accordingly into blocks, as in Eq. Assume 

PNS-in(O) = "P N s-inPs(0)'P NS _ in 7^ °- Then the factor 
7Yns is called a decoherence-free (or noiseless) subsystem 
if the following condition holds: 

Tr in { PNS _ in (i)} = U NS Tr in {p NS _ in (0)}U^ s , (49) 

where Uns is a unitary matrix acting on 7Yns ■ 

Definition 6 Perfect initialization (DF subsystems): 
p' = and p out — in Eq. 

Definition 7 Imperfect initialization (DF subsystems): 
p' and/or p out in Eq. are non-vanishing. 

According to Definition |5j a quantum state encoded 
into the Hns factor at some time t is unitarily related to 
the t = state. The factor 7ii n is unimportant, and hence 
is traced over. Clearly, a NS reduces to a DF subspace 
when 7ij n is one-dimensional, i.e., when H- m = C. 

We now present the necessary and sufficient conditions 
for a NS and later we show that the algebra-dependent 
definition, Eq. Qllfl. is a special case of this generalized 
form. In stating constraints on the form of the Kraus 
operators, below, it is understood that in addition they 
must satisfy the sum rule J2 a ElE a = I, which we do 
not specify explicitly. 
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Theorem 4 Assume imperfect initialization. Then a 
subsystem H-^s i n the decomposition 7ig = 7Yns ® W n © 
7i out is decoherence-free ( or noiseless ) with respect to CP 
maps iff the Kraus operators have the matrix representa- 
tion 



E, 



U® C Q 
B„ 



Corollary 2 Assume perfect initialization. 
Kraus operators have the relaxed form 



(50) 



Then the 



E„ 



B„ 



(51) 



We note that this result has been recently derived 
from an operator quantum error correction perspective 
in Ref. 20] . Note again that there is a trade-off between 
the quality of preparation and the amount of leakage that 
can be tolerated, a fact that was not noted previously 
for subsystems, and has important experimental impli- 
cations. 

As discussed above, the original definition of a NS 
was based on representation theory of the error alge- 
bra. Here we have argued in favor of a more comprehen- 
sive definition, based on the quantum channel picture. 
Let us now state explicitly why our result is more gen- 
eral. Indeed, in the algebraic approach one arrives at 
the representation l|13|l of the Kraus operators, namely 
E« = ©jgjlnj ® G a: j. However, it is clear from 
Eq. ||inj that our channel-based approach leads to a form 
for the Kraus operators that includes this latter form as 
a special case, since it allows for the off-diagonal block 
A Q . The representation (|13[l of the Kraus operators does 
agree with Eq. Q5U|). but in that case we do not need to 
assume initialization inside the NS, so that again, our 
result is more general than the algebraic one. 



2. Markovian Dynamics 

As in the CP-map based definition of a NS, we need 
to trace out the 7i ln factor, here in order to obtain the 
dynamical equation for the subsystem factor: 



dt 



dt 



= Tr in { 



d^NS-inAS^NS-in 



dt 



} 



Tr in {7> N s_ in (-|[H S! ps] + ^2F Q p 5 Ft 



h " ,r "' ' 2 
-Ft F aPs - PS FiF a )r^_. m }. 



(52) 



Definition 8 The factor 7Yns is called a decoherence- 
free ( or noiseless ) subsystem under Markovian dynamics 
if a state subject to Eq. \5ty) . undergoes continuous uni- 
tary evolution: 



p m = i[M,p NS ] 



(53) 



where M is Hermitian. 

Clearly, again, a NS reduces to a DF subspace when 
Win is one-dimensional, i.e., when H m = C. 

Our goal is to find necessary and sufficient conditions 
such that Eq. (|52|l leads to Eq. I|53|) . In the case of per- 
fect initialization, since it does not involve Ti. out , Eq. I|52|) 
is meaningful only if the system remains in the subspace 
Hns ® Win- An analysis of Eq. l(52*j l reveals that this 
leakage-prevention goal is achieved by imposing the con- 
straints stated in the following theorem, proven in Ap- 
pendix |X2j 



Theorem 5 Assume perfect initialization. Then a sub- 
system Hns in the decomposition Hs — Ti-NS'^'HinO'Hout 
is decoherence-free (or noiseless) with respect to Marko- 
vian dynamics iff the Lindblad operators have the matrix 
representation 



F — 



Ins ® C a A Q 







(54) 



and the system Hamiltonian (including a possible Lamb 
shift) has the matrix representation 



H NS (8) Iin+lNS ® H in H 2 

Hs ! Hi 



L 2 ±1 3 

where H m is constant along its diagonal, and where 



(55) 



H 2 



■^(l NS ®Ct)A 



(56) 



Eqs. I)55 [l . I)56 |l are new additional constraints on the 
Lindblad operators (compared to Ref. [10() which must 
be satisfied in order to find a NS. 

If, on the other hand, we allow for imperfect initializa- 
tion, we find a different set of conditions: 

Theorem 6 Assume imperfect initialization. Then a 
subsystem 7Yns in the decomposition Hs = Hns <8> Win © 
7i ou t is decoherence-free ( or noiseless ) with respect to 
Markovian dynamics iff the Lindblad operators have the 
matrix representation 



B„ 



(57) 



and the system Hamiltonian (including a possible Lamb 
shift) has the matrix representation 



H 



H 



NS 



Iin+lNS ® H in 







H, 



(58) 



IV. PERFORMANCE OF QUANTUM 
ALGORITHMS OVER IMPERFECTLY 
INITIALIZED DFSs 

In this section we discuss applications of our gen- 
eralized formulation of DFSs to quantum algorithms. 
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As mentioned above, a major obstacle to exploiting 
decoherence-free methods is the unrealistic assumption 
of perfect initialization inside a DFS. Removing this con- 
straint enables us to perform algorithms without perfect 
initialization, while not suffering from information loss. 
We separate the role of an initialization error in the al- 
gorithm (i.e., starting from an imperfect input state), 
from the effect of noise in the output due to environment- 
induced decoherence. Thus we hrst quantify an error en- 
tirely due to incorrect initialization (Ai ca k below), then 
compare the DFS situations prior and post this work, by 
relating them to Ai ca k- 

1) Initialization error in the absence of decoherence: 
Assume no decoherence at all, that the initial state is 



actual / 



Pi 92 
P\ P3 



p — (0)= 

while the ideal input state is fully in the DFS: 



P o 




(59) 



(60) 



Further assume that the algorithm is implemented via 
unitary transformations U = Udfs © Idfs- 1 - > applied to 
T^dfs- I n general this will lead to an output error in the 
algorithm, which can be quantified as 



A loak = ||Up actual (0)U t -Up idcal (0)U t || 
Udfs(pi - p)Up, FS U DF sp 2 



p\v 



DFS 



P3 



(61) 



where ||-|| denotes an appropriate operator norm. This 
error appears not because of decoherence but because 
of an erroneous initial state. This is a generic situation 
in quantum algorithms, which is not special to the DFS 
case: Eq. (|59fl is generic in the sense that one can view 
the DFS block as the computational subspace, with the 
other blocks representing additional levels (e.g., a qubit 
which is embedded in a larger Hilbert space). Methods 
for correcting such deviations from the ideal result exist 
(leakage elimination H3,|5l|), but are beyond the scope 
of this paper. 

2) Initialization error in the presence of decoherence: 
Assume that the input state is imperfectly initialized, as 
in Eq. (|59|l . and in addition there is decoherence, i.e., 



P 



actual 



(62) 



with the Kraus operators given by Eq. I|19f) [the form 
compatible with decoherence-free evolution starting from 
p actual (0)]. Prior to our work it was believed that for an 
imperfect initial state of the form p actual (0), leakage due 
to the components pi and p^ would cause non-unitary 
evolution of the DFS component. Thus instead of an 
error U D fs(/Oi - /o)Ud FS in the DFS block of Eq. @TJ, it 
was believed that one had £{p\) — UdfspUd FS where £ 
is an appropriate superoperator component. This would 



have led to a reduced algorithmic fidelity, A| cak < Ai ea k- 
However, we now know that even for an initial state of 
the form p actual (0), when the Kraus operators are given 
by Eq. H19fl the actual algorithmic fidelity is still given 
by Aieak, since in fact the evolution of the DFS block is 
still unitary. 

The above arguments apply when imperfect initializa- 
tion is unavoidable but one knows the component p\. 
A worse (though perhaps more typical) scenario is one 
where not only is imperfect initialization unavoidable, 
but one does not even know the component p\. In this 
case the above arguments apply in the context of al- 
gorithms that allow arbitrary input states. Almost all 
the important examples of quantum algorithms are now 
known to have a flexibility of this type: Grover's algo- 
rithm |22j was the first to be gene ralized to allow for 
arbitrary input states, first pure [2J, |23, |26j , then mixed 
[27j : Shor's algorithm [2lJ can run efficiently with a sin- 
gle pure qubit and all other qubits in an arbitrary mixed 
state 2§|; a similar result applies to a class of inter- 
esting physics problems, such as finding the spectrum 
of a Hamiltonian |2^; the Deutsch-Josza j2^| algorithm 
was generalized to allow for arbitrary input states |30j| . 
and a similar result holds for an algorithm that performs 
the functional phase rotation (a generalized form of the 
conventional conditional phase transform) [3lJ . Most re- 
cently it was shown that Simon's problem and the period- 
finding problem can be solved quantumly without initial- 
izing the auxiliary qubits |32^ | . 

For algorithms that do not allow arbitrary input states, 
one could still make use of the flexibility we have intro- 
duced into DFS state initialization, provided it is possible 
to apply post-selection: one modifies the output error of 
algorithm by observing whether the measurement out- 
come came from the DFS block or not (this could be 
done, e.g., via frequency-selective measurements, simi- 
lar to the cycling transition method used in trapped-ion 
quantum computing |59^. 



V. DECOHERENCE FREE SUBSPACES AND 
SUBSYSTEMS IN NON-MARKOVIAN 
DYNAMICS 

A. Decoherence Free Subspaces 

In Ref. [3i| a new class of non-Markovian master equa- 
tions was introduced. The following equation was derived 
as an analytically solvable example of this class: 



^ = -i[U s , p s ] +cj* dt'k{t') exp(jM)p s (t - *') 



(63) 



where C is Lindblad super-operator and k(t) represents 
the memory effects of the bath. The Markovian limit is 
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clearly recovered when k(t) oc S(t). 2 

Some examples of physical systems which can be de- 
scribed by this master equation are (i) a two-level atom 
coupled to a single cavity mode, wherein the memory 
function is exponentially decaying, k(t) — e~ xt [^3, and 
(ii) a single qubit subject to telegraph noise in the par- 
ticular case that ||£|| <C 1/i, whence Eq. <|H>5|l reduces to 
Ps — C J Q * dt'k(t')p(t — t') [6(j. It is interesting to inves- 
tigate the conditions for a DFS in the case of dynamics 
governed by Eq. (|(j3fl . and to compare the results with 
the Markovian limit, k(t) oc S(t). We defer proofs to 
Appendix IA 31 and here present only the DFS-condition, 
stated in the following theorem (note that, similarly to 
the Markovian case, we consider here a continuous-time 
DFS). 



Theorem 7 Assume imperfect initialization. Then a 
subspace 7Ydfs is decoherence free iff the system Hamil- 
tonian Hg and Lindblad operators F a have the matrix 
representation 



He 



Hdfs 








H 



DFS- 1 



c Q I 
B a 



(64) 



These conditions are identical to those we found in the 
case of Markovian dynamics with imperfect initialization 
- cf. Theorem [5] This fact provides evidence for the 
robustness of decoherence-free states against variations 
in the nature of the decoherence process. 

Interestingly, the conditions under the assumption of 
perfect initialization differ somewhat when comparing 
the Markovian and non-Markovian cases: 



Corollary 3 Assume perfect initialization. Then a sub- 
space 7Ydfs is decoherence free iff the system Hamilto- 
nian H5 and Lindblad operators F Q have the matrix rep- 
resentation 



H s = 
F Q = 



H 



DFS 

H 



DFS J 



c a T A c 
B n 



and c* A Q = 0. 



(65) 
(66) 



B. Decoherence Free Subsystems 

We now consider the NS case. The dynamics governing 
a NS is derived by tracing out 7Yi n : 



dp 



NS 



Ot 



dTr in {p s } _ dps_ 

dt ~ ini dt 1 
TVin-J— i[H s ,ps] 

+£ f dt'k(t')exp(Ct')p s (t-t')} (67) 
Jo 



Theorem 8 Assume imperfect initialization. Then a 
subsystem Hns in the decomposition H.s = H-ns ® Ti-in © 
Hout is decoherence-free (or noiseless) with respect to 
non-Markovian dynamics [Eq. i)6'3)) / iff the Lindblad op- 
erators and the system Hamiltonian have the matrix rep- 
resentation 



I NS ® C a 
B Q 

H NS ® Iin+lNS O H in 

H„„ 



(68) 
(69) 



Note that this form is, once again, identical to the 
Markovian case with imperfect initialization (cf. Theo- 
rem EJ) . 

However, as in the DFS case, the conditions are slightly 
different between Markovian and non-Markovian dynam- 
ics if we demand perfect initialization: 

Corollary 4 Assume perfect initialization. Then a sub- 
system H^s in the decomposition Hs — T^NS^^inffi'^out 
is decoherence-free (or noiseless) with respect to non- 
Markovian dynamics [Eq. iff the Lindblad operators 
and the system Hamiltonian have the matrix representa- 
tion 



Ins ® C a A a 
B Q 

5](Ins®C!)A q = 0, 

rf _ , H NS <g) Iin+lNS ® Hj, 





H ou t 



VI. SUMMARY AND CONCLUSIONS 



(70) 
(71) 

(72) 



Compared to the Markovian case (Theorem^ , the dif- 
ference is that now the off-diagonal blocks of the Hamilto- 
nian must vanish, whereas in the Markovian case we had 
the constraint [Eq. PdfsH s P dfs = -i^ a c*A Q . 



2 We note that Ref. I.'j.'jj contains a small error: the Markovian limit 
is recovered for k(t) = S(t) only if the lower limit in Eq. 1631 is— t. 
This change can easily be applied to the derivation of Ref. I.'i.'l . 



We have revisited the concepts of decoherence-free sub- 
spaces and (noiseless) subsystems (DFSs), and intro- 
duced definitions of DFSs that generalize previous work. 
We have analyzed the conditions for the existence of 
DFSs in the case of CP maps, Markovian dynamics, and 
(for the first time) non-Markovian continuous-time dy- 
namics. Our main finding implies significantly relaxed 
demands on the preparation of decoherence-free states: 
the initial state can be arbitrarily noisy. If, on the other 
hand, the initial state is perfectly prepared, then almost 
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arbitrary leakage from outside the DFS into the DFS can 
be tolerated. 

In the case of Markovian dynamics, if one demands 
perfect initialization, our findings are of an opposite na- 
ture: we have shown that then an additional constraint 
must be imposed on the system Hamiltonian, which im- 
plies more stringent conditions for the possibility of ma- 
nipulating a DFS than previously believed. We have pre- 
sented an example to illustrate this fact. 

We have also shown that the notion of noiseless sub- 
systems, as originally developed using an algebraic ap- 
proach, admits a generalization when it is instead devel- 
oped from a quantum channel approach. 

Our results have implications for experimental work 
on DFSs. and in particular on quantum algorithms over 
DFSs tLfltLfl- It is now known that a large class of quan- 
tum algorithms can tolerate almost arbitrary preparation 
errors and stillprovide an advantage over their classical 
counterparts [MlllllllllllllMlMllllll. The 
relaxed preparation conditions for DFSs presented here 
are naturally compatible with this approach to quantum 
computation in noisy systems. This should provide fur- 
ther impetus for the experimental exploration of quan- 
tum computation over DFSs. 
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as Hns <8 H in © Hout- From the NS definition, Eq. (@HJ), 
we have 

Tr in {U ® I(P NS -inPs(0)^ s _ in )Ut ® 1} = 

Tr in {^ (7> NS -inE Q ) Ps(0)(Et ^ S _.J}. (A2 ) 

a 

Let us represent the Kraus operators in the same block- 
structure matrix- form as that of the system state, i.e., 
corresponding to the decomposition Hs = Hns <S> Hi n © 
^out, where the blocks correspond to the subspaces 
Hns © Hin (upper-left block) and H ou t (lower-right 
block). Then 



PS 
E„ 



Pi P2 

pi pa 

P A 

Da B Q 



(A3) 
(A4) 



with appropriate normalization constraints, considered 
below. Equation (|A2|) simplifies in this matrix form as 



Tr in {U ® I Pl Ut <g> 1} = Tr in {J] P aj0l pt 

a 

+V aP2 A-t + A aJ 4P+ + A aP3 Al}, (A5) 

which must hold for arbitrary ps(0). To derive con- 
straints on the various terms we therefore consider special 
cases, which yield necessary conditions. First, consider 
an initial state ps(0) such that pi = 0. Then, as the LHS 
of Eq. (|A5|) is independent from p 3 , the last term must 
vanish: 



(A6) 



APPENDIX A: PROOFS OF THEOREMS AND 
COROLLARIES 

Here we present proofs of all our results above. We 
shorten the calculations by starting from the NS case 
and obtain the DFS conditions as a special case. 



1. CP Maps 

a. Arbitrary Initial State 

Assume the system evolution due to its interaction 
with a bath is described by a CP map with Kraus op- 
erators {E Q }: 



ps(*)=^E Q p s (0)E. 



(Al) 



Note that here ps is an operator on the entire system 
Hilbert space Tls, which we assume to be decomposable 



Further assume pi = <8> |«')(«'|. Note that the par- 
tial matrix element (j'|P Q |i') is an operator on the 7Yns 
factor, Then Eq. (|A5II reduces to 

\i){i\ = ]T [U+OlPalO] [<i'|P£li'>U] • (A7) 

a,j' 

Taking matrix elements with respect to a state or- 
thogonal to \i), yields: 



o = ^|(^|[ut(/|p Q |0]K)| : 

X | [Ut<j'|P a |i'>] \i)=0, 



a 



(A8) 



which, in turn implies that [U^(j'|P a |z')] \i) is propor- 
tional to |z), i.e., 



[<j'|Pa|i'>]|i>cxU|i>. 



(A9) 



Since are arbitrary this condition implies that 

the submatrix P Q must be of the form P a = U ® C Q . 
Substituting P a — U <£> C a into Eq. (|A5f) we have 
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Tr in {U ® IpxUt ® 1} = Tr in {J2 a U g> C Qj0l Ut g> Ct }, 
so that 



When (0) = 7- , d/ 3 s(0)Pd we have that /? 2 = and 
p 3 = and Eq. I|A5|I reduces to 



Tr in {pi} =TV in {^I NS ®C Q piI NS ®Ct}. (A10) Tr in {U <g> I^U* $ 1} = Tr in {^ P Q piP^}. (A16) 



Now supp ose pi = Y^iji'j' hji'j'\i)(j\®> then from 

Eq. ijA"T0)l we find 

iji' 

E (k'\C a \i')(j'\cUk').(All) 

iji' j'k' a. 

Using J2k' \ k ')( k '\ = !in, Eq. jMTt becomes 

E A «^"N)0'l = E ^|i)(j|(i'|EctC a |0. 

zjV iji' j' a 



It follows that 



E ^"a^a — li) 



(A12) 



(A13) 



Next consider the normalization constraint 
J2 a EqE q = I for the Kraus operators, together with 
the additional constraints we have derived (A a = 0, 
P a = V®C a ): 

EPiP a + D^D Q = I NS ®I in 

a 

i NS <g> e c l c « + E D « D « = Ins 

a a 

(A14) 

But, from Eq. (|A13|) we have J] Q P a P a = r NS ® lin- 
Therefore D a = 0. 

Taking all these conditions together finalizes the ma- 
trix representation of the Kraus operators as 



E Q — 



U«C a 
B„ 



(A15) 



For a scalar C a we recover the DFS condition 1)190. 
These considerations establish the necessity of the rep- 
resentation (|A15|I : it is simple to show that this repre- 
sentation is also sufficient, by substitution and checking 
that the NS and DFS conditions are satisfied. Therefore 
we have proved Theorems ^ an d El 



b. Perfect Initialization 

We now prove Corollaries and [21 for DF-initialized 
states of the form ps{0) = 'PdPs{0)'Pd- Thus, we have to 
prove that D a = in Eq. (|A4|) . 



The argument leading to the vanishing of the A a 
[Eq. HA6fl ] then does not apply, and indeed the A a 
need not vanish. However, the arguments leading to 
P a = U <g> C Q and £ a P^Pa = Ins ® lin do apply. 
Hence D„ = 0. 



2. Markovian Dynamics 

a. Arbitrary Initial State 
Consider Markovian dynamics 



dps 
dt 



-i[U s ,p s ] + J2 F a PsFi 

a 

—Ft F aPs - ipsFt F Q , (A17) 

with the following matrix representation of the various 
operators: 



, Pi Pi 

PS = f 

Pi P3 



He 



Hi H 2 \ ( P Q A Q 

Ht H 3 1 



Da B ( , 



(A18) 



Then we find the dynamics of the NS block to be 

dpNS _ dTrjnjpi} _ 

dt dt 
-*'Tr in {[H llPl ]} - iTr in {(H 2 4 - p 2 H 2 )} + 

Tr in {E PaPiPj* + A a p|pt + P a p 2 Al + A a p 3 Al 

a 

~\ E( P « P « + D ^ D «Vi + ( P l A - + DtB Q )4 

E Pi( p l p « + D ^ D «) + P2(At P« + BlD a )} 

(A19) 

The right-hand side of this equation must be independent 
of p2 and /?3, for any matrices p 2 and p%. Therefore the 
term Aq^A^ implies A Q = 0. Collecting the remaining 
terms acting on p\ from the left yields Trj n {(— zH 2 — 
P>iB a )pl} = 0. Together we have 



A a = 0, iH 2 + E D « B a = 0. 



(A20) 
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This reduces Eq. ljA19|) to 

dpm #Tr in {pi} 



dt dt 

-iTtin[Hi, Pi] + Tr in PcPiPl 

a 

4 Tr -E{( P « P « + D i D «)"°i} (A21) 



Consider the initial state p\ = pns <8 with € 



dpNS 

at 



-i[<*'|Hi|i'),A>Ns] 



■ £<j'|p Q |i'>p NS <i'|pt|i') 

a 

|Ews,((*'iPii/)(j'iPoN') 

a 

(z'|Di|/)(/|D Q |0)} 



(A22) 



Let /5ns = IV*) I with ijj arbitrary and apply (ift ), 
such that = 0, to Eq. I|A22|I . denoting P Q v y = 



El(^|Pa,^}| 2 = o. 



(A23) 



Since this identity must hold for all ip and ■0- L , we find 
that P a y,j' — c a,i',j' Ins , which implies that P a = 
Ins <8 Cfjj. Moreover, by definition of a NS, there exists a 
Hermitian matrix Hns such that pns obeys a Schrodingcr 
equation, dp^s/dt = — i[Hns,Pns]- Therefore the non- 
Hermitian term DjjD Q in Eq. (|A21|) must vanish, 
implying that D Q = 0. 

Combining these results with Eq. (|A20|I yields 



ffftin{pi} 

dt 



= -iTr ln {[H 1)Pl ]} 
= -i[H NS ,p NS ] 



(A24) 



This identity can be realized iff Hi = Hns ® Iin+lNS <8> 
Hi n . Therefore the NS conditions are obtained as 



H 

¥ 



Hns ® Iin+lNS ® H in 
H 3 

Ins®C? 
B„ 



(A25) 



The DFS condition is a special case of l|A22|l . with 
dim(Win) = 1. This concludes the proof of Theorems [3 
and|Hl 



6. Perfect Initialization 
Now consider perfect initialization: 
Ps -- 



Pi =0 
= =0 



(A26) 



This is just the case of an arbitrary initial state con- 
sidered above, with p 2 — and p 3 = in Eq. I|A19|I . 
This then yields the dynamics of pns as being given by 
Eq. (|A2ip . Repeating the derivation following Eq. (|A21|) 
we conclude again that D Q = 0, P Q = Ins ® Cf n and 
Hi = H NS ® Iin+ lNS ® H in . 

Note that Eq. (|A20J| now docs not apply (it was ob- 
tained assuming nonzero p 2 , P3), i.e., we cannot conclude 
that A Q and H2 vanish. This implies that that dps/dt 
has a non-zero off-diagonal elements, which, using the 
master equation (|A17|I . we calculate to be: 

upper right block: 

i Pl H 2 + £ P Q piDt - ipi(PlA Q + D+ B ) 

= i Pl H 2 - i pi 53(I NS 8) C^ n + )A Q 

a 

bottom right block: VjD a piD^ = 0. 

a 

To prevent the appearance of corresponding off-diagonal 
blocks in pg, we must therefore demand 



H 2 + ^(Ins ® C^)A a = 0, 



(A27) 



which is Eq. 1)56(1 . The DFS case is obtained with 
dim(?ii n ) = 1. This concludes the proof of Theorems |3| 
and 



3. Non-Markovian Dynamics 

The derivation of the conditions for decoherence- 
freeness in the case of non-Markovian dynamics is some- 
what different from the other two cases we have consid- 
ered, because of the appearance of the nonlocal-in-time 
integral in the master equation: 

^ = -i[K s , p s ] +£ f dt'k(t') exp(£t')ps(i - 

(A28) 

In order to find necessary conditions on the structure of 
H5 and C consider the case of small t, expand 



Ps 



(*) =J2 tn Ps ^(0), k(t) = E * m fc (m) (0), (A29) 



n=0 



m=0 



and substitute into Eq. (|A28|) . The constant (t°) term 
yields 



pW(0) = -i[H SlPS (0)]. 



(A30) 



The terms involving t 1 yield, after Taylor-expanding 
exp(£i'): 

2p£°(0) = -m s ,p^(0)]+k(0)Cp s (0). (A31) 
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Thus the solution of Eq. (|A28|I up to first and second 
order in time is: 

Ps(t) = p s (0)-it[H s ,p s (0)}+O(t 2 ), (A32) 
ps(t) = ps{0)-it[H s ,p s (0)} 
t 2 



-{-[H s ,[H 5 ,p 5 



■k(0)Cp s (0)} + O(t 3 ). 

(A33) 



a. Arbitrary Initial State 

Consider once again the matrix representations as in 
Eq. (|A18|) . Substituting these expressions into the first 
order equation (|A32(I , the pi it) block yields 

Pns(*) = «is(0)-*m lll {[Hi > pi(0)]} 

-ziTr in {H 2 p 2 (0) - p 2 (0)H 2 } => 
H 2 = 0, Hi =H NS ®I in + I NS ®H in . 

(A34) 

Continuing to second order, Eq. I|A33(I . the NS block is 
found to be 

pNs(t) = Pns(O) - «i[H NS , pns(O)] 



Hns^ fH 



NS 



p NS (0)]]+Tr in {2fc(0)^P aPl pt 



-fc(0) ^(PlP Q + Dt D Q )p! + (PlA a + Dt B a )pl 

a 

-k(0) Pi(PlPa + Dt D Q ) + pa (At P Q + BlD Q )}. 

a 

(A35) 

The first three terms correspond to unitary evolution, but 
the remaining terms are essentially identical to the case 
of Markovian dynamics and must be made to vanish, just 
as in Eq. (|A19|) . The same arguments used there apply 
and consequently 



InsOC" 
B„ 



(A36) 



The conditions (|A34|) . i|A36(l are necessary and sufficient 
for unitary evolution of the NS block under our non- 
Markovian master equation. The DFS case is obtained 



with dim(7ii n ) = 1. This concludes the proof of Theo- 
rems [7| and 00 



b. Perfect Initialization 



Assume 



Ps(0) 



(A37) 



then from the first order equation (|A32|) , the NS block is 
found to satisfy 



pNs(i) = PNs(0)-ifIhn{[Hi,p(0)]}=> 

Hi = H N s <8 Iin + Ins ® H in . (A38) 

To second order in time [Eq. l|A33|l ]: 

PNs(t) = Pns(0) - it[H NS , pns(O)] 
-^[H ns ,[H ns ,pns(0)]] 

+ ^Tr in {-H 2 H 2 p(0) - p(0)H 2 H 2 
+2k(Q) P aP Pl - (Pt P„ + DlD Q )p(0) 

OL 

-p(0)(ptp Q + DtD a )}, 

(A39) 

which is again similar to the Markovian case. Similar 
logic therefore yields H 2 = D Q = 0, and hence 



Ins ® C a A D 
B n 



(A40) 



Here we should notice that the density matrix ps(0) 
has an off-diagonal element p(0) ^ a (Pj,A Q + Dj,B Q ) = 
p(0)^ Q PtA a . This term must vanish, for otherwise 
Ps(t) has non-zero off-diagonal elements. Summarizing, 
we have 



Ins ® C Q A a 




H 



B°)' E( In s® C «) A «=° 

/ a 
H NS ® Iin+lNS ® H in 







H„ 



(A41) 



The DFS case is obtained with dim(7ij n ) = 1. This con- 
cludes the proof of Corollaries and 0] 



[1] R. Landauer, Proc. Roy. Soc. London Ser. A 353, 367 

(1995) . 

[2] W.G. Unruh, Phys. Rev. A 51, 992 (1995). 

[3] S. Haroche and J.M. Raimond, Physics Today 49, 51 

(1996) . 

[4] L.-M Duan and G.-C. Guo, Phys. Rev. A 57, 737 (1998). 



[5] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 
(1997). 

[6] DA. Lidar, I.L. Chuang and K.B. Whaley, Phys. Rev. 

Lett. 81, 2594 (1998). 
[7] DA. Lidar, D. Bacon, J. Kempe, and K.B. Whaley, Phys. 

Rev. A 63, 022306 (2001). 



15 



[8] E. Knill, R. Laflamme and L. Viola, Phys. Rev. Lett. 84, [32 
2525 (2000). [33 
[9] S. De Filippo, Phys. Rev. A 62, 052307 (2000). 
[10] J. Kempe, D. Bacon, D.A. Lidar, and K.B. Whaley, Phys. [34 

Rev. A 63, 042307 (2001). 
[11] C.-P. Yang and J. Gea-Banacloche, Phys. Rev. A 63, [35 
022311 (2001). 

[12] P.G. Kwiat, A.J. Berglund, J.B. Altepeter, and A.G. [36 

White, Science 290, 498 (2000). 
[13] D. Kielpinski, V. Meyer, M.A. Rowe, C.A. Sackett, W.M. [37 

Itano, C. Monroe, and D.J. Wineland, Science 291, 1013 

(2001) . [38 
[14] E.M. Fortunato, L. Viola, J. Hodges, G. Teklemariam, 

and D.G. Cory, New J. Phys. 4, 5 (2002). [39 

[15] L. Viola, E.M. Fortunato, M.A. Pravia, E. Knill, R. 

Laflamme, D.G. Cory, Science 293, 2059 (2001). [40 

[16] M. Mohseni, J.S. Lundeen, K.J. Resch, A.M. Steinberg, [41 
Phys. Rev. Lett. 91, 187903 (2003). [42 

[17] J. Ollerenshaw, D.A. Lidar, and L.E. Kay, Phys. Rev. [43 
Lett. 91, 217904 (2003). 

[18] D.A. Lidar, K.B. Whaley, in Irreversible Quantum Dy- 
namics, Vol. 622 of Lecture Notes in Physics, edited by [44 
F. Benatti and R. Floreanini (Springer, Berlin, 2003), p. 
83, eprint quant-ph/0301032 [45 

[19] J. Preskill, in Introduction to Quantum Computation and 

Information, edited by H.K. Lo, S. Popescu and T.P. [46 
Spiller (World Scientific, Singapore, 1999). [47 

[20] DW. Kribs, R. Laflamme, D. Poulin, and M. Lesosky, 
(2005), eprint quant-ph/0504189 

[21] PW. Shor, SIAM J. on Comp. 26, 1484 (1997). [ 

[22] L.K. Grover, Proceedings of the 28th Annual ACM Sym- [49 
posium on the Theory of Computing (ACM, New York, 
NY, 1996), p. 212. [50 

[23] D. Deutsch and R. Jozsa, Proc. Roy. Soc. London Ser. A 

439, 553 (1992). [51 

[24] D. Biron, O. Biham, E. Biham, M. Grassl and D.A. Lidar, [52 
in Quantum Computing & Quantum Communications; [53 
First NASA International Conference; selected papers, [54 
QCQC'98, Vol. 1509 of Lecture Notes in Computer Sci- 
ence, edited by C. P. Williams (Springer, Palm Springs, [55 
1998), pp. 140-147, eprint quant-ph/9801066 [56 

[25] E. Biham, O. Biham, D. Biron, M. Grassl and D.A. Lidar, 

Phys. Rev. A 60, 2742 (1999). [57 

[26] E. Biham, O. Biham, D. Biron, M. Grassl, D.A. Lidar, 

and D. Shapira, Phys. Rev. A 63, 012310 (2001). [58 

[27] E. Biham, D. Kenigsberg, Phys. Rev. A 66, 062301 

(2002) . [59 
[28] S. Parker and M.B. Plenio, Phys. Rev. Lett. 85, 3049 

(2000). 

[29] E. Knill and R. Laflamme, Phys. Rev. Lett. 81, 5672 

(1998). [60] 
[30] D.P. Chi, J. Kim, and S. Lee, J. Phys. A 34, 5251 (2001). 
[31] J. Kim, S. Lee and D.P. Chi, J. Phys. A 35, 6911 (2002). 



D.P. Chi, J.S. Kim, and S. Lee, eprint quant-ph/0504173 

A. Shabani and D.A. Lidar, Phys. Rev. A 71, 020101(R) 
(2005). 

P. Zanardi and M. Rasetti, Mod. Phys. Lett. B 11, 1085 
(1997). 

D. A. Lidar, D. Bacon and K.B. Whaley, Phys. Rev. Lett. 
82, 4556 (1999). 

G. M. Palma, K.-A. Suominen and A.K. Ekert, Proc. Roy. 
Soc. London Ser. A 452, 567 (1996). 

L.-M Duan and G.-C. Guo, Phys. Rev. Lett. 79, 1953 

(1997) . 

P. Zanardi and F. Rossi, Phys. Rev. Lett. 81, 4752 

(1998) . 

P. Zanardi and D.A. Lidar, Phys. Rev. A 70, 012315 
(2004). 

L. Viola, E. Knill, Phys. Rev. A 68, 032311 (2003). 
S. Nakajima, Prog. Theor. Phys. 20, 948 (1958). 
R. Zwanzig, J. Chem. Phys. 33, 1338 (1960). 

H. -P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, Oxford, 
2002). 

K. Kraus, States, Effects and Operations, Fundamental 
Notions of Quantum Theory (Academic, Berlin, 1983). 
V. Gorini, A. Kossakowski and E.C.G Sudarshan, J. 
Math. Phys. 17, 821 (1976). 

G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 
R. Alicki and K. Lendi, Quantum Dynamical Semigroups 
and Applications, No. 286 in Lecture Notes in Physics 
(Springer- Verlag, Berlin, 1987). 

E. B. Davies, Rep. Math. Phys. 11, 169 (1977). 

D.A. Lidar, A. Shabani, and R. Alicki, Chem. Phys., in 
press (2005), eprint quant-ph/0411119 
D.A. Lidar, Z. Bihary, and K.B. Whaley, Chem. Phys. 
268, 35 (2001). 

P. Zanardi, Phys. Rev. A 57, 3276 (1998). 

N.P. Landsman, eprint math-ph/9807030 

DW. Kribs, Proc. Edinburgh Math. Soc. 46, 421 (2003). 

J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 

68, 580 (1992). 

N. Gisin, I.C. Percival, J. Phys. A 25, 5677 (1992). 

M. Plenio and P. Knight, Rev. Mod. Phys. 70, 101 

(1998). 

L.-A. Wu, M.S. Byrd, D.A. Lidar, Phys. Rev. Lett. 89, 
127901 (2002). 

M.S. Byrd, D.A. Lidar, L.-A. Wu, and P. Zanardi, Phys. 
Rev. A 71, 052301 (2005). 

D. J. Wineland, C. Monroe, W. M. Itano, D. Leibfried, 

B. E. King, and D. M. Meekhof, J. of Res. of the Na- 
tional Inst, of Standards and Technology 103, 259 (1998), 
http://nvl.nist.gOv/pub/nistpubs/jres/103/3/cntl03-3.htm 
S. Daffer, K. Wodkiewicz, J.D. Cresser, J.K. Mclver, 
Phys. Rev. A 70, 010304 (2004). 



